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Abstract. The radial scale length of the exponential component of the disc of the Milky Way has been determined 
in the near infrared. We have used the TMGS (Two Micron Galactic Survey) database which contains positions and 
K-magnitudes of about 700.000 stars distributed in serveral regions along the galactic plane. From those we have 
selected areas more than 5° off the plane to minimize the effect of extinction and the contributions of the young 
disc; and with longitudes ranging from 30° to 70°, to avoid contaminations from the central bulge, the bar and the 
molecular ring, in the inner end, and from the local arm, the warp and the truncation of the disc, in the outer end. 
We observed stars with magnitude 9 < tuk < 10. The use of the NIR K-band also reduces the effect of extinction. In 
the observation region, mjr = 9.5 mag stars are K2-K5 III stars, with an absolute magnitude that is nearly constant, 
which also greatly simplifies the problem. We have obtained the value 2.1 ± 0.3 kpc for the radial scale length, which 
is a typical value when compared with external galaxies of similar types. 
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1. Introduction 

It is well known that normal discs of spiral galaxies have a component which can be fitted to an exponential (Freeman, 
1970) with a radial scale length h. There are two basic reasons to undertake a new determination of the parameter h 
in the Milky Way: 

a) There is a very large scatter in the values found in the literature. For instance, Kent, Dame and Fanzio (1991) 
after reviewing previously reported values, found that they range from 1 to 5.5 kpc, and proposed an intermediate 
value of 3 kpc. van der Kruit (1986) obtained 5.5 kpc; Lewis and Freeman (1989) gave a value of 2.5 kpc, and so 

' did Robin, Creze and Mohan (1992), among others. This uncertainty is also too large in the NIR. Mikami and Ishida 
(1981) obtained a value of 1 kpc, Jones et al. (1981) 2 kpc, and Eaton et al. (1984) 3 kpc. Values of h in the NIR seem 
to be lower, but even this fact remains uncertain, despite its interest in the dynamic evolution of the disc. 

b) Even if h is simply a fitting parameter for a distribution with frequent deviations from an exponential, it 
has became a natural radial length unit, very useful when comparing spiral galaxies with different sizes or unknown 
distances, and physically more significant than true length of angle units. The position of any morphological feature 
or typical lengths of disc phenomena are usually and hopefully expressed taking h as unit. Some of these features in 
our own galaxy may be considered normal or exceptional depending on the value of this normalizing length. 

Probably the large scatter of data is a result of managing a large amount of data concerning different regions of 
the sky, which are difficult to analyse because of the contribution of many other complex components, such as bulge, 
star formation ring associated with the gas ring, bar, spiral arms, warp and so on. However, a proper choice of the 
wavelength and the region of the sky may render this problem a very simple one, which implies a more confident 
obtention of the parameter h 

This search specifically deals with, and aims at an improved determination of the radial scale length of the Milky 
Way disc. 

We have chosen the K-band, in which extinction is very low. This also implies that the observed distribution closely 
resembles the true one. In addition, we have taken measurements at b = 5°, rendering extinction negligible. We then 
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surveyed a galactic longitude range from 30° to 70°. This practically eliminates the contribution of bulge, ring and bar 
(Hammersley et al., 1994) which take place for I < 30°, and the contribution of the Local Arm, warp and truncation 
(Porcel et ah, 1997), which take place for I > 70°. The contribution of the other spiral arms (not the Local Arm) is 
minimized taking 6 = 5°, as they remain closer to the plane. Another advantage of working in the K-band is that the 
observed stars in the 9-10 magnitude range are practically only K2-K5 III stars, through the surveyed sky zone. This 
fact was realised after taking into account the previous analyses of Ruelas-Mayorga (1991), Wainscoat et al. (1992), 
Calbet et al. (1995) and our own initial calculations (Porcel, 1997). The light coming from these stars is 80 % of the 
total light from all stars in this magnitude range. This greatly simplifies the problem as the luminosity function can 
be approximated by a Dirac's delta function. 

Observations were carried out with the Carlos Sanchez 1.5 m telescope at the Teide Observatory, as part of the 
Two Micron Galactic Survey (TMGS) project (Garzon et ah, 1993; Hammersley et ah, 1994; Calbet et ah, 1995). 
Interpolation techniques were used to obtain a regular mesh from the series of discrete scans (Porcel, 1997; Porcel et 
al. 1997). 



2. Results 

As explained above we consider K2-K5 III stars. The luminosity function 4> is defined as usual, i. e. <p(M)dM gives 
the number of stars per pc 3 with an absolute K magnitude in the interval [M, M + dM] . It is here approximated 
by a Dirac's function, i. e. 4>{M) — S(M — Mo) centred at Mq = — 2.5mag. Its true dispersion is only 0.6 mag. The 
effect of a non- vanishing dispersion on absolute magnitude could be important, as the resulting value for the scale 
length could be affected by some sort of Malmquist bias. We have performed several numerical calculations using a 
gaussian distribution for the luminosity function, obtaining similar results. The final value of the estimated error is not 
affected either, as it mainly arises from the fitting process. Even though these kind of effects are difficult to evaluate, 
the approximation of the luminosity function by a S — function was found to be a satisfactory assumption. In the 
fundamental equation of star counts (e. g. Mihalas and Binney, 1981; Gilmore and Reid, 1983; Calbet et ah, 1995) the 
density function is also taken into account. It was assumed that: 

R ~ R Q |z| 

n(R,z) = e s ~ (1) 

where n(R, z) is the density function, i. c. the number density of stars taking its local value as unit, R is the galac- 
tocentric distance, and z is the vertical coordinate. We take Rq — 8.5 kpc. The vertical exponential profile and the 
value of h z = 200 pc were adopted from Wainscoat et al. (1992). We have repeated these calculations for different 
values of h z and we have obtained similar results, i. e. neither the proposed value of h nor its estimated error are 
affected. Actually h is found to be noticeably independent of h z . We therefore think that this choice for h z cannot be 
an important source of errors. 

Defining as usual A(l,b), as the number of stars per squared degree with magnitude between 9 and 10, the 
fundamental equation of star counts becomes 

r 10 r°° r- r \z\ 

A = lu dm CS(M-M )e s ~r 2 dr (2) 



where u> is the solid angle, C is a constant and r is the distance from the Sun. This expression is easily transformed 
into 

A(l) = C 3 e — (3) 

where 

#0 = (R 2 Q + r 2 Q cos 2 b - 2R r o cos b cos I) 1/2 (4) 
C3 is another constant, and 

m — Mn + 5 

r = 10 s = 2.5 kpc (5) 

As b is taken to be constant, Ro is a function only of I, and therefore A is a function only of /. Equation (3) predicts 
a linear function relation between log .4 and Rq- The results are plotted in Figure 1. In this figure, error bars represent 
Poisson errors. After fitting (using a standard chi-squared method as described in Press et ah, 1992), we obtained 

h= (2.1 ±0.3) kpc (6) 
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Fig. 1. Log of the number of K=9-10 magnitude stars per squared degree as a function of Ro(l) defined in the text. 
3. Conclusions 

The value we propose, h = (2.1 ± 0.3) kpc has been obtained by choosing the most appropriate conditions for this 
task. This value is similar to typical values for other galaxies (de Grijs and van der Kruit, 1996; Peletier et al., 1994; de 
Jong, 1996). This value seems quite low when compared with the obtained values in the optical (see for instance van 
der Kruit, 1986; and Barteldress and Dettmar, 1994) but it is well known than the exponential scale lengths decrease 
with increasing wavelength (de Grijs and van der Kruit, 1996; Peletier et al., 1994; de Jong 1996, Tully et al., 1996). 
Peletier et al. (1994) found that the ratio hs/hx is in the range 1.2-2 in a sample of 37 Sb's and Sc's edge-on galaxies, 
with this ratio increasing with axis ratio, showing that it is an effect that is mostly due to extinction. In this case 
our determination of the radial scale length of the Milky Way is compatible with values in optical bands as large as 
4 kpc. Our value is virtually free of extinction effects, so it is a proper determination of the true scale length of the 
mass distribution of the Milky Way disc. A similar value (2.3 ± 0.1 kpc )has been obtained by Ruphy et al. (1996). 
Fux and Martinet used indirect methods based on the asymmetric drift equation, and also obtained a similar value 
(2.5lg;g kpc). We therefore conclude that with respect to the radial scale length, our galaxy is a typical one. 
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